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Abstract 
Introduction. Properties of filled composites based on polytetrafluoroethylene allow them to work at high contact 


pressures, reciprocating nature of shear loads, and in a wide temperature range. Due to this, they are used as antifriction 
layers of bearing parts with ball segment. To simulate the mechanical behavior of such materials under operating 
conditions, adequate constitutive equations of elasto-viscoplasticity and methods of their identification according to the 
data of basic experiments are needed. 

Materials and Methods. The tensor-linear model of elasto-viscoplasticity should be identified according to the data of 
tests on free compression of samples. They were subjected to loading up to a maximum deformation of 10 %, allowed 
to remain, unloaded, and then, a similar loading cycle up to 160 MPa under constrained compression was carried out. 
The experiment with a composite based on polytetrafluoroethylene filled with 40 wt. % fine bronze, was conducted at 


room temperature. Tests for constrained compression were performed for two values of the strain rate, and for free 
compression — for three values of the strain rate in the range of 10°—10~s''. For the description, two models of 


elasto-viscoplasticity were considered, representing modifications of Swain and Kletschkowski’s models and 
corresponding to the connection of a viscoelastic or elastic nonlinear viscous element with a plastic or endochronic 
element. An integral operator with a Kohlrausch kernel was considered as a viscoelastic element. 

Results. The results of the constrained compression tests made it possible to separate the elastic relationship of 
volumetric deformations and average stresses from the constitutive relations. The data of free compression cycles at 
different strain rates were used to determine the material constants of the model. For this purpose, an efficient search 
algorithm based on the simplex method of minimizing the discrepancy was implemented. Both models discovered the 
importance of the plastic component (independent of the deformation rate) for a qualitative description of the stress 
cycling that accompanied the cyclic deformation, as well as their dependence on the strain rate. 

Discussion and Conclusions. Both models of elasto-viscoplasticity described correctly the behavior of the studied 
fluorocomposite under loading conditions close to the operating conditions of the antifriction layers of the bearing parts 
with ball segment. They can be considered as a basis for their further generalization, taking into account the dependence 


on temperature. 


Keywords: elasto-viscoplasticity, constitutive equations, filled fluorocomposite, identification, free and constrained 


compression, trapezoidal loading. 


Acknowledgments. The authors are grateful to A. V. Khokhlov, leading researcher of Institute of Mechanics, 


Lomonosov Moscow State University, for a fruitful discussion of the work. 


©D. S. Petukhov, A. A. Adamov, I. E. Keller, 2022 


Petukhov D. S., et al. Selection and Identification of a Model of Elasto-Viscoplasticity of the Filled Fluorocomposite according to Free 


Funding information. The research is done with the financial support from RFFI and Perm Krai within the frame of 
research project no. 20-48-596012. 


For citation. D.S.Petukhov, A.A. Adamov, I. E. Keller. Selection and Identification of a Model of Elasto- 
Viscoplasticity of the Filled Fluorocomposite according to Free and Constrained Compression Tests. Advanced 
Engineering Research, 2022, vol. 22, no. 3, pp. 180—192. https://doi.org/10.23947/2687-1653-2022-22-3-180-192 


Introduction. Properties of filled composites based on polytetrafluoroethylene (PTFE) allow them to work in a 
wide temperature range at high contact pressures and reciprocal sliding on the counterbody. Due to this, they are used as 
antifriction layers of bearing parts with a ball segment'*. PTFE as a base for antifriction polymer composites in a wide 
temperature range shows high stability of properties and exceptionally low coefficients of friction, sliding and rest. The 
introduction of fillers helps to significantly increase the wear resistance, stiffness and yield strength, as well as reduce 
the creep of the composite without marked increase in the friction coefficients. The rheological properties of the 
antifriction material should provide feasibility during the operation of the product in the temperature range from less 


than —50 to +50 °C: 
—at calculated values of the compressive load of at least 60 MPa; 
—at its peak values up to 150 MPa; 
—at a normalized range of cyclic shear deformations. 


The predictability of changes in the properties and thickness of the antifriction layer under its hinge mounting also 


matters. 


During the operation of the materials under consideration, it is required to assess the performance of antifriction 
layers of the support parts with a ball segment, to predict the resource for the entire service life (up to 50 years). To do 
this, you need to construct and identify adequate constitutive equations of elasto-viscoplasticity according to the data of 
basic experiments. Models describing the rheological behavior of filled fluorocomposites are also important for 
calculating the operational properties of shaft seals and rods in high-power diesel engines*. To calculate the behavior of 
products under load during installation (for short periods of time), it is important that the models correctly describe the 
plastic (elastoplastic, viscoplastic) properties of the material. To calculate the behavior of products during operation (for 
long periods of time), it is required to adequately represent the phenomena of relaxation, creep, and cyclic creep 


(ratcheting). 


Rheological models of filled fluorocomposites describing cycles of loading, holding, and unloading in wide ranges 
of time and strain rates are rarely found in the literature. We can specify only work [1] that meets all these requirements. 
It proposes a model corresponding to the parallel connection of endochronic and nonlinear viscous elements. This 


model describes all uniaxial tensile tests* at room temperature. It is about: 


— loading curves up to 0.08 strain with strain rates 10“ —107s"'; 


— relaxation curves for 10 h with constant strain of 0.03—0.08; 


— hysteresis loops in the load cycle up to 0.1 strain, short-term hold and unloading with subsequent hold. 


Mechanics 


' EAD 050009-00-0301. Spherical and cylindrical bearing with special sliding material made of fluoropolymer. EOTA. Product regulation. URL: 
https://www.nlfnorm.cz/en/ehn/6189 
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One-dimensional constitutive equations containing seven material constants were written, as well as their tensor 
generalization for the geometrically linear case. The authors point out that the endochronic plasticity model makes it 
possible to describe a continuous change in plastic deformation during loading and unloading. This is typical for 
polymer materials and is not typical for metals. Note that in [1], the endochronicity parameter was not used, which 
allowed regulating the characteristic range of transition to the plastic state. In addition, the paper’ describes similar 
curves for temperatures up to 120°C, and strain rates in the range of 10*-100 s", relaxation curves with a sharp 


increase in temperature up to 100 °C, and relaxation curves under uniaxial compression. 


A tensor geometrically nonlinear model is presented in [2]. It describes the behavior of filled fluorocomposites 
under various monotonic and cyclic loading histories in various stress states and at arbitrary temperatures. A 
viscoplastic rheological element in a series connection with a viscoelastic one is taken into account. On the whole, this 
model satisfactorily describes the curves corresponding to the experiments [1], but it is much more complicated and 
requires fourteen material constants. In [3], the model [4] is used to describe microindentation cycles of a filled 
fluorocomposite. The rheological block diagram is similar to [2], but it takes into account the relaxation time spectrum 
of the viscoelastic element. In [3], the model [4] is used to describe microindentation cycles of a filled fluorocomposite. 
The rheological block diagram is similar to [2], but it takes into account the relaxation time spectrum of the viscoelastic 
element. 

The model [5] was designed to describe features of the elastoplastic behavior of solid polymers under cyclic uniaxial 
loading. The evolutionary equations for its internal variable, stresses and plastic deformations, had nonlinear cross- 
connections, therefore, the model did not correspond to any rheological structural scheme. The measure of internal time 
characteristic of the endochronic theory of plasticity was used. Twenty-one constants were required to identify the 
model. 

The nonlinear one-dimensional Maxwell-type viscoelastic plasticity model [6-9] contains only two material 
functions. It describes incomplete recovery of deformations in the load-unloading cycle and cyclic creep (ratcheting) 
with asymmetric loading cycles. Nonlinear viscoelasticity models with fractional integro-differentiation operators are 
characterized by flexibility with a small number of material constants and functions [7]. It is necessary to investigate the 
applicability of such nonlinear models to describe the behavior of polymer materials in the loading and unloading cycle 
with holds [10-13]. Papers [14, 15] are devoted to the identification of linear models of viscoelasticity based on the 
results of indentation tests. In [16, 17], the dynamic mechanical analysis results are used for this purpose, but the 
significant nonlinearity of the materials studied in this work does not allow using this technique. Ratcheting of PTFE 


composites is experimentally investigated in [18, 19], but there is no mathematical description of the results. 


Paper [8] presents a method for identifying a model of isotropic elastoplastic behavior of filled composites based on 
PTFE and ultrahigh molecular weight polyethylene based on experimental data on free and constrained compression. 
Due to the peculiarities of the rheology of plastics [9], the methods [10] are applicable to them. This approach makes it 


possible to determine: 
— the elastic volume compressibility function; 
— the hardening function for shear elastoplastic properties within the framework of the theory of plastic flow. 


Thus, the operating conditions of the products in question, the mechanical properties of polymer materials, and the 


presented literature analysis serve as the basis for further research. 


Let us highlight the basic test for free compression of cylindrical samples of 20x20 mm. It consists of load cycles up 
to 0.1 deformation, a short-term hold of fifteen minutes, and unloading with a subsequent hold of fifteen minutes. In this 


way, data can be determined to describe elastic, plastic and viscous components of deformation, relaxation and reverse 


> Material and Computational Analysis of PTFE Seals. 
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creep with residual deformations. A similar definition of hardness is provided by the standard®. Loading is performed at 
three strain rates in the range of 10°-10°s". Besides, the samples are tested for constrained compression to a stress of 
160 MPa at two strain rates from the same range. All tests are repeated for values from the operating temperature range 
of the polymer material. These data should be sufficient to identify a suitable model of elasto-viscoplastic compressible 


material with nonlinear elastic compressibility. 


First of all, we are interested in the possibilities of describing experimental data at room temperature with the 


simplest elasto-viscoplasticity models of the type proposed in [1] and [4], which differ in structural schemes. 


Materials and Methods. A composite based on PTFE PN 90 filled with 40 wt. % of fine bronze was chosen as the 
simulated material. Cylindrical samples with a diameter and height of 20 mm were made by pressing and sintering a 


mixture of powders. The following experiments were conducted to select and then identify the model: 
— for constrained compression with a single nonzero component of the strain tensor; 


— for free compression with a single nonzero component of the stress tensor. 


In the first case, friction along the lateral boundary was excluded by lubrication, in the second — by a thin film of 
PTFE and grease. 


2,000 4,000 6,000 ts 


a) 
o, MPa 


150 


100 


50 


b) 


Fig. 1. Cyclic tight compression test: a — program; b — hysteresis curve 
The experiments were aimed at identifying irreversible volumetric deformations and volume-elastic properties of the 
material. The first test was cyclic constrained compression with increasing amplitude (Fig. 1 a). It set the value of the 
deformation rate €=2.5x10~°s”" and the pause between unloading and loading for 300 s. During the first few cycles, a 


residual volume strain equal to approximately 5x10~° accumulated (Fig. 1). It changed slightly when the sample was held in 


the unloaded state. 


® GOST 4670-2015 (ISO 2039-1:2001). Plastics. Determination of hardness. Ball indentation method. Interstate Council for Standardization, 
Metrology and Certification. Moscow: Standartinform; 2016. 10 p. (In Russ.) 
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Sg and 


Further, compression tests were carried out up to 160 MPa with a significantly slower strain rate € =2.5x10 


subsequent unloading (Fig. 2 a): 
— of undeformed sample (1); 
— of the same sample repeatedly (2); 


—of the sample after the cyclic test shown in Figure 1 (3). 


o, MPa 
150 Test no.1 
— Test no.2 
— Test no.3 
100 
50 
0 
0.01 0.02 0.03 0.04 & 
— Testno.1 
— Test no.2 
— Test no.3 


0.001 0.002 0.003 0.004 € 
b) 


Fig. 2. Slow tight compression tests: a — loading diagrams; b — initial and final sections of diagrams 


Tests of the same sample were carried out in a day. 

As can be seen in Figure 2 b, in the experiment with a new sample, the maximum residual strain of about 3.5x10° is 
observed. On repeated testing, it amounted to 1x10”. In the third test, there is no residual strain. The properties of the 
material under cyclic constrained compression are accompanied by a transient irreversible process and are stabilized in 
several cycles. 

The bulk properties were determined from the third test (Fig. 2 a), in which there are no residual strains. In this 
paper, bulk strains are assumed to be nonlinearly elastic. The data is described by a quadratic dependence: 

o=M(sje; M(e)=a,+a,6; a, =2.4hPa; a, =22.1 hPa. (1) 

Here, 6,6 — axial components of the Cauchy stress tensors and logarithmic strains; M — constrained compression 


modulus, which is related to the bulk compression modulus K_ by the ratio: 


K(e,)(3K(e,) +E 
Mee, )=3 ReIGKE) +B.) (2) 
9K(e,)—-E,, 
where E,, — equilibrium Young's modulus (independent of bulk strain); €, =€;, € — components of the logarithmic 


strain tensor. 
Value E,, =690 MPa was obtained with the slope of the curve o(€) under uniaxial compression in the vicinity of 


o = 1MPa and relaxation 8.65 %. 
Using (1), (2) and 9.38 % correction for the compliance of the machine in the tight compression test, it is possible to 


obtain the desired dependence K(é,) , which differs slightly (maximum by 0.03 %) from linear: 
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0,, =K(e, Jey; K(e,)=b,+b,; by, =2.1hPa; b, =22.2 hPa, (3) 


where 6,, =+6,, 6; — components of the Cauchy stress tensor. 


i> OF 

Experiments to determine the relationship between average stress and volume strain were carried out on samples 
previously tested by constrained compression with irreversible volume strain. The same applies to the free compression 
experiments described below and required to determine the relationship of the deviatory parts of stress and strain 
tensors. As a result, the model reflected the properties of the compacted material. In production, processes associated 
with irreversible bulk strain occur under the crimping of a layer of antifriction polymer material during the assembly of 
the structure of the support part and do not affect the further behavior of the material. 

Uniaxial compression tests were carried out to determine the relationships between the deviatory parts of stress and 


strain tensors. Figure3 shows the results of experiments for three values of strain rates 
&=4.63x10°s', 2.35x10°s', 2.27x10°s''. Further, the data from these experiments will be interpreted in the 


framework of a one-dimensional model. When formulating a three-dimensional model, the elastic volume part can be 
excluded from these relations. 


o, MPa 


Experiments at different speeds 
15 — ¢=4.63x10° 
— €=2.35x10" 
— €=2.27x10" 
10 
5 
0 ZZ | 


0.02 0.04 0.06 008 010 « 
a) 


0.10 Experiments at different speeds 
, — £=4.63x10° 
0.08 - €=2.35x10° 
— €=2.27x10° 
0.06 
0.04 
0.02 


5 10 15 20. 25 1x10°, s 


Experiments at different speeds 
15 - ¢=4.63x10° 
- €=2.35x10° 
— €=2.27x10% 
10 
5 
0 


5 10 15 20 «25 1x10°, s 
c) 


Fig. 3. Uniaxial compression experiment with three different speeds: 
a — loading diagrams; b — deformation histories; c — stress histories 


Mechanics 


185 


http://vestnik-donstu.ru 


186 


Advanced Engineering Research 2022. V. 22, no. 3. P. 180-192. ISSN 2687-1653 


When choosing a suitable model, it turned out that the experimental data in Figure 3 were not described by the 
relations of a linear viscoelastic medium: 


o(t)=[R(t-r)e(r)d7 = Ry [ec -{r¢ pena} (4) 


where o(t), €(t) — stress and strain histories used to describe elastomers and polymers [10]. 
There have been attempts to use expressions for relaxation functions R(t) or relaxation rates 7“(t) with a different 


number of parameters: the sum of exponents, the Robotnov fractional-exponential function’, the Koltunov kernel’, the 
kernel from the work of T. L. Smith’, the Kohlrausch kernel [10]: 


—-rR iS ie Thee _ Aexp(-Bt”). 
RO=R, +). R, exp(—a,t); M(t) = At LT (ain) T(t) = : s 
R,-R.,, 7 : 
RD Rots Gaye KO=Re Re Rae tay"). 


Here, R,,, Ry, R,;, @, 8, T — material constants; [ — gamma function, as well as monotonic approximation by 
cubic splines at eight points, whose coordinates served as parameters. 
Models were considered, in which, in addition to viscous ones, plastic irreversible deformations were taken into 


account. Two basic models of this type were investigated. 
Basic model | (Fig. 4a) is a combination of models [4] and [10] and assumes a sequential connection of three 


elements. 


Fig. 4. Block diagrams of basic models: a — model 1; b — model 2 


1. Linear viscoelastic element with Kohlrausch kernel [10]: 


: a 
€,(t) = J T(t—r)6(r)dr, ‘7 (t) = 7 (z em Jos ‘) (6) 
where R, R,,, 7), @ — material constants, the first two have the physical meaning of instantaneous and equilibrium 
modules. 
2. Plastic element with linear hardening: 
. k6(t), C26, Vvo<S-o, 
s10=| 0, -o,<o<o, ° 7) 
where o,,, k — material constants. 
3. Linear-viscous element: 
é,()=1 'o(t), (8) 
where 1 — material constant. 
The final expression for strains: 
e(t) =€,(t) +€,(t)+e,(t). (9) 


7 Rabotnov YuN. Polzuchest' ehlementov konstruktsii. Moscow: Nauka; 1966. 752 p. (In Russ.) 

* Koltunov MA. Singulyarnye funktsii vliyaniya v analize relaksatsionnykh protsessov. Prochnost' i plastichnost'. Moscow: Nauka; 1971. P. 640-645. 
(In Russ). 

° Smith TL. Ehmpiricheskie uravneniya dlya vyazkouprugikh kharakteristik i vychisleniya relaksatsionnykh spektrov. Vyazkouprugaya relaksatsiya v 
polimerakh. Moscow: Mir; 1974. P. 44—56. (In Russ). 


Petukhov D. S., et al. Selection and Identification of a Model of Elasto-Viscoplasticity of the Filled Fluorocomposite according to Free 


Model (6)-(9) is solved with respect to strains, therefore, it is convenient to have stress histories for its 
identification o(t) (Fig. 3 c). 

The block diagram of basic model 2 (Fig. 4 b), borrowed from [1], is a parallel connection of two elements. 

1. Endochronic plastic element with nonlinear elastic part: 


é()= SO lac) ; 6,(t) = Asign(e(t) —e,(t))In(1 +Ble(t)—e, (t) ), (10) 
where €, — an internal parameter that has the meaning of irreversible strain of the element, Y,A,B — material 
constants. 

2. A pair of sequentially connected nonlinear-viscous (2) and linear-elastic (3) elements: 
k 
a,(t)= #sin(o, (0) ; (1) 
7] 0 
0,(t) =0,(t) =ce,(t), (12) 
where sign — argument sign; n, 0,, k , c — material constants. 
The resulting expression for stresses: 
o(t) =0,(t)+o,(t). (13) 


Model (10)—(13) is solved with respect to stresses, therefore, it is convenient to have strain histories for its 
identification g(t) (Fig. 3 c). Expression o,(t) in (13) is solved as the solution to a system of algebraic differential 
equations (10), and o,(t) — as a system solution (11), (12) taking into account g(t) =¢,(t)+¢,(t). 

Research Results. Residual minimization between experimental strain histories (Fig. 3 b) and model predictions 
allowed us to find seven material constants of base model 1: 

R =895 MPa;R,, = 205 MPa; t, = 1.66x10's; a = 0.64; 
6, =10.2 MPa; k = 0.0095 MPa; 1 =9.28x10°s. 


The simplex method is used as a search procedure. Comparison of the calculated and experimental data is shown 


(14) 


in Figure 5. 
& 
Slow experiment 
0.10 — Experiment 
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0.08 e Residual points 
0.06 
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Quick experiment 
0.10 =o — Experiment 
— Model 


0.08 e Residual points 


250 500 750 1,000 1,250 1,500 t,s 
c) 
Fig. 5. Comparison of experimental data and base model 1 predictions: 
a— &=4.63x10°s';b— €=2.35x10°s";c— €=2.27x10" s* 


Residual minimization between the experimental stress histories (Fig. 3 b) and the model predictions allowed us to 


detect seven material constants of base model 2: 


Y = 20.3MPa; A = 6.01 MPa; B= 61.8; 


15 
1 =6.81-10° s; 6, = 0.66 MPa; k = 1.08; s = 681 MPa. (15) 


Here, strain histories were considered to be given (Fig. 3 b). Comparison of the calculated and experimental data is 


shown in Figure 6. 


o, MPa 
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o, MPa 


Quick experiment 


— Experiment 
— Model 
e Residual points 


250 500 750 1,000 1,250 1,500 tx10°, 


c) 
Fig. 6. Comparison of experimental data and predictions of base model 2: 


a— €=4.63x10°s";b— €=2.35x10°s';c— €=2.27x10" Ss" 
Discussion and Conclusions. The considered one-dimensional models were generalized to a spatial formulation 
that assumed isotropy of properties and immutability of the main axes of the strain tensor under loading. To do this, the 


Cauchy stress tensors and logarithmic deformations were decomposed into spherical and deviatoric parts: 
1 
o=o0,I+58, € mae 


The relationship between the ball parts 6,, and €, is given by expression (3). The relationship between the 


m 


deviatoric parts s and e is given by the relations generalizing (6)—(9) or (10)-(13) depending on the selection of the 
base model. These expressions should be rewritten in terms of the histories of the deviatoric parts s(t) and e(t), 


excluding the ball part from uniaxial compression. 

AS a Spatial generalization of the relations (7) for the plastic element of model 1, the plastic flow law associated with 
von Mises yield criterion, with isotropic linear strain hardening was considered. The connection of deviators will look 
like this: 

e(t) =e,(t) +e, (t) +e, (t); 


e,(t) = | 7(t-1)8(r)dr; M)== (i 1 Jos( (t/t,)*); 


t 


e, =2k ay e,(t) = fé, (od; Gy = fee, 26s; 
e, 0 


-l 
e,(t)=n s(t), 
where colon means the operation of convolution of tensors A: B= A,B, given by the components in an orthonormal 


basis. 
The relations of the endochronic plasticity theory (10) of model 2 are also generalized in the spirit of von Mises 
plasticity, which provides the following entry for deviators: 


S=S,+8,;e=e, +e; 


3 
6, = 56S; &=,/26:6; 5 =dev{ Y'éinn,} é, = Asign (A, )In(1+B 
i=l 


di 


; 


k 


31 fo) S 
2 2 ae =. | 8a. eae & 
é, 5 sh 3 0, =/38, :8); 
n 99 /| S 
=n iy 


where “dev” means a tensor deviator, A,, n; — eigenvalues and eigenvectors of the tensor e—e, =, n, nj. 
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It should be noted that the authors [1] used a model option with an endochronic parameter [14] equal to one. It 
seems appropriate to consider this parameter as a material constant, since its influence on the description of transients 
after changes of loading modes is known. 

In the framework of this work, the experiments were carried out on constrained and free compression of a composite 
based on polytetrafluoroethylene filled with finely dispersed bronze with a mass fraction of 40 %. At the same time, 
loading, unloading and holding were carried out with different strain rates in the range of 10°-10™ s"'. The experiments 
on constrained compression were carried out at pressures corresponding to the peak loads of the antifriction layer in 
hinges with a spherical segment. The inability of the standard linear viscoelastic model to describe test data for several 
strain rates was found. Two families of seven-constant models with plastic and nonlinear viscous structural elements 
that could describe the data of basic experiments were selected and identified. The model with an endochronic element 
showed great accuracy in describing the data at a high strain rate. The model with a plastic element can be improved if 
the nonlinearity in the viscous element and nonlinear strain hardening are taken into account, which will require an 
increase in the number of material constants. The data of these models were generalized to the spatial case for an 


isotropic material and strain histories with a constant orientation of the trihedron of the main axes. 


References 
1. Kletschkowski T, Schomburg U, Bertram A. Endochronic Viscoplastic Material Models for Filled PTFE. 
Mechanics of Materials. 2002;34:795—808. https://doi.org/10.1016/S0167-6636(02)00197-7 


2. Bergstrém JS, Hilbert Jr. LB. A Constitutive Model for Predicting the Large Deformation Thermomechanical 


Behavior of Fluoropolymers. Mechanics of Materials. 2005;37:899—9 13. https://doi.org/10.1016/j.mechmat.2004.09.002 


3. Stan F, Munteanu AV, Fetecau C. Viscoelastic Characterization of Polytetrafluoroethylene (PTFE) Polymer by 


Sharp Indentation. AIP Conference Proceedings. 201 1;1315:221—226. https://doi.org/10.1063/1.3552445 


4. Stan F, FetecauC. Study of Stress Relaxation in Polytetrafluoroethylene Composites by Cylindrical 


Macroindentation. Composites Part B: Engineering. 2013;47:298—307. 10.1016/j.compositesb.2012.11.008 


5. Mencik J, Li Hong He, Swain MV. Determination of Viscoelastic-Plastic Material Parameters of Biomaterials by 
Instrumental Indentation. Journal of the Mechanical Behavior of Biomedical Materials. 2009;2:318—325. 


https://doi.org/10.1016/j.jmbbm.2008.09.002 


6. Drozdov AD, J de C Christiansen. Cyclic Elastoplasticity of Solid Polymers. Computational Materials Science. 


2008;42:27-35. https://doi.org/10.1016/j.commatsci.2007.06.002 


7. Khokhlov AV. Comparative Analysis of Creep Curves Properties Generated by Linear and Nonlinear Heredity 
Theories under Miulti-Step Loadings. Mathematical Physics and Computer Simulation. 2018;21:27-51. 
https://doi.org/10.15688/mpcm.jvolsu.2018.2.3 


8. Khokhlov AV. Applicability Indicators and Identification Techniques for a Nonlinear Maxwell-Type Elasto- 
Viscoplastic Model using Multi-Step Creep Curves. Herald of the Bauman Moscow State Technical University. Natural 
Sciences. 2018;6:92—112. http://dx.doi.org/10.18698/1812-3368-2018-6-92-112 

9. Khokhlov AV. Identification Techniques for the Nonlinear Maxwell-Type Viscoelastoplastic Model Using Creep 
Recovery Curves. Problems of Strength and Plasticity. 2018;80:238—254. https://doi.org/10.32326/18 14-9 146-2018-80- 


2-238-254 


Petukhov D. S., et al. Selection and Identification of a Model of Elasto-Viscoplasticity of the Filled Fluorocomposite according to Free 


10. Ogorodnikov EN, Radchenko VP, Ungarova LG. Mathematical Models of Nonlinear Viscoelasticity with 
Operators of Fractional _Integro-Differentiation. PNPRU ~~ Mechanics Bulletin. 2018;2:147-161. 
https://doi.org/10.15593/perm.mech/2018.2.13 


11. Adamov AA. Experimental Verification and Identification of the Models of Isotropic Body Showing Elastic 
Volume Compressibility of Disperse-Filled Composites Based on Polytetrafluorethylene and Ultra-High-Molecular- 
Weight Polyethylene. Konstruktsii iz kompozitsionnykh materialov. 2013;2:28—37. 

12. Goldman AYa. Ob"emnoe deformirovanie plastmass. Leningrad: Mashinostroenie; 1984. 232 p. (In Russ.) 

13. Adamov AA, Matveenko VP, Trufanov NA, et al. Metody prikladnoi vyazkouprugosti. Ekaterinburg: URO 
RAN; 2003. 411 p. (in Russ.) 

14. Cheng L, Xia X, Yu W, et al. Flat-Punch Indentation of Viscoelastic Material. Journal of Polymer Science. Part 
B: Polymer Physics. 2000;38: 10-22. 10.1002/(SICD 1099-0488(20000101)38:1<10::AID-POLB2>3.0.CO;2-6 


15. Lu H, Wang, B, Ma J, etal. Measurement of Creep Compliance of Solid Polymers by Nanoindentation. 


Mechanics of Time-Dependent Materials. 2003;7:189—207. https://doi.org/10.1023/B:MTDM.0000007217.07156.9b 
16. Toth LF, P De Baets, Szebényi G. Thermal, Viscoelastic, Mechanical and Wear Behaviour of Nanoparticle 


Filled Polytetrafluoroethylene: A Comparison. Polymers. 2020;12:1940. https://doi.org/10.3390/polym1209 1940 


17. Gavrilenko SL, Shilko SV. Attestatsiya lineinoi vyazkouprugoi modeli antifriktsionnogo polimernogo 
kompozita po rezul'tatam uskorennykh ispytanii na relaksatsiyu. Teoreticheskaya i prikladnaya mekhanika. 
2017;32:155—158. (In Russ.) 

18. Konova EM, Ostrer SG, Khatipov SA. Vliyanie prirody napolnitelya na fiziko-mekhanicheskie svoistva 
radiatsionnykh modifikatsii kompozitov na osnove politetraftorehtilena. Plasticheskie Massy. 2011;5:40—43. (In Russ.) 

19. Wenjuan Xu, Hong Gao, LiLanGao, et al. Tensile Ratcheting Behaviors of Bronze Powder Filled 
Polytetrafluoroethylene. Frontiers of Chemical Science and Engineering. 2013;7:103—109. 


https://doi.org/10.1007/s11705-013-1315-8 


20. Kadashevich Yu, Pomytkin S. Endochronic Model of Plasticity Generalizing Sanders's Theory. Scientific 
Letters of Rzeszow University of Technology. Mechanics. 2014;31:539-547. https://doi.org/10.7862/RM.2014.57 


Received 08.08.2022 
Revised 29.08.2022 
Accepted 31.08.2022 


About the Authors: 
Petukhov, Dmitriy S., lead engineer of the Laboratory of Nonlinear Mechanics of Solids, Institute of Continuous 


Media Mechanics UrB RAS (1, Academika Koroleva St., Perm, 614018, RF), ResearcherID, ScopusID, ORCID, 


petuhovds @ mail.ru 


Adamoy, Anatoliy A., leading researcher of the Laboratory of Nonlinear Mechanics of Solids, Institute of 
Continuous Media Mechanics UrB RAS (1, Academika Koroleva St., Perm, 614018, RF), Dr.Sci. (Phys.-Math.), 
ScopusID, ORCID, adamov @icmm.ru 


Mechanics 


191 


Advanced Engineering Research 2022. V. 22, no. 3. P. 180-192. ISSN 2687-1653 


Keller, Ilya E., head of the Laboratory of Nonlinear Mechanics of Solids, Institute of Continuous Media Mechanics 


UrB RAS (1, Academika Koroleva St., Perm, 614018, RF), Dr.Sci. (Phys.-Math.), ResearcherID, ScopusID, ORCID, 


kie @icmm.ru 


Claimed contributorship: 


D. S. Petukhov: modification of models and identification of their material constants according to experimental data. 


A. A. Adamov: development of methodology and conducting experiments. I. E. Keller: literature review and model 


selection. 


Conflict of interest statement 


The authors do not have any conflict of interest. 


All authors have read and approved the final manuscript. 


http://vestnik-donstu.ru 


192 


